Initial state preparation with dynamically generated system-environment correlations 
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The dependence of the dynamics of open quantum systems upon initial correlations between the 
system and environment is an utterly important yet poorly understood subject. For technical con- 
venience most prior studies assume factorizable initial states where the system and its environments 
are uncorrelated, but these conditions are not very realistic and give rise to peculiar behaviors. One 
distinct feature is the rapid build up or a sudden jolt of physical quantities immediately after the 
system is brought in contact with its environments. The ultimate cause of this is an initial imbalance 
between system-environment correlations and coupling. In this note we demonstrate explicitly how 
to avoid these unphysical behaviors by proper adjustments of correlations and/or the coupling, for 
setups of both theoretical and experimental interest. We provide simple analytical results in terms 
of quantities that appear in linear (as opposed to afhne) master equations derived for factorized 
initial states. 



I. OPEN SYSTEM INITIAL CORRELATIONS 

An open quantum system is a quantum system 'S' that 
interacts with some environment 'E' whose degrees of 
freedom have been coarse grained (colloquially, 'traced 
out' or 'integrated over' ). The unitary evolution of the 
combined system -f environment, 'C (for combined or 
closed), is generated by the Hamiltonian 

He = Hs + He + Hi , (1) 

where Hi denotes the system-environment interaction. 
Specifying the initial of the combined system is neces- 
sary to determine the open-system dynamics. The most 
common choice is to assume factorized initial states for 
the system and environment^ 

Pc((I) = Ps(0)®Pe(0), (2) 

where Ze{13) denotes the partition function of the free 
(noninteracting) environment and T = 1//3 is the tem- 
perature of the environment, which acts here as a thermal 
reservoir. 

When considering environments with a large number 
of high-frequency modes and characterized by a UV fre- 
quency cutoff A, such a factorized initial state (chosen 
for mathematical simplicity) unfortunately engenders un- 
physical behavior such as a sudden jolt in physical quanti- 
ties near the initial time (this was analyzed in some detail 
in Ref . [ ] ) or spurious cutoff sensitivity of certain system 



Even after resolving issues of renormalization, the system will 
typically be displaced by a finite amount of the order of the 
induced damping within a very short time of the order of the 
inverse UV cutoff of the environment. 



correlators (see Ref. p] and references therein) . This kind 
of initial conditions assumes that an uncorrelated sys- 
tem and environment are instantaneously coupled with 
non-vanishing strength. The pathological behavior arises 
because the factorized initial state contains a number of 
highly excited energy states of the full Hamiltonian (in- 
cluding the interacting Hamiltonian) , even when the ini- 
tial reduced states of the system and environment are not 
highly excited in the free theory, and it is a reflection of 
the high-frequency modes of the environment quickly be- 
coming correlated with the system within a time of order 
1/A. 

The next most common choice of initial state (see 
Ref. [ ]) has been to consider system deformations or 
measurements of the global equilibrium state of the com- 
bined system 'C, with density matrix 

Pc(0) = X^O:,^e-^H-O„, (4) 

where the O and O' operators are restricted to act on 
the system. However, this method still gives rise to jolts 
for sufficiently general deformations or measurements [ ] , 
which can be understood as a consequence of altering the 
state of the system instantaneously [■")] . 

To cure or tame these drastic effects, especially in 
the context of linear systems, the following procedure 
has been suggested: a) force the system by a constant 
amount, b) wait for it to relax into the displaced equi- 
librium state, and then c) release the force [■>]. Alterna- 
tively and in order to generate interesting coherent su- 
perposition states for the system, one can start with the 
equilibrium state of the combined system and act on the 
system, but for a non- vanishing time [•)]. Essentially we 
view the problem as an imbalance between initial corre- 
lations and initial coupling strength; the imbalance can 
be countered on either side. We also believe that the 
most natural resolution should be a dynamical prepara- 
tion which relies upon equilibration [ !, followed by an 
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additional preparation of the system for a finite time ["] . 
Our key contribution is showing that this can be achieved 
while still taking advantage of the simpler analytical re- 
sults obtained when deriving the master equation for a 
factorizcd initial state, without the need to introduce in- 
homogeneous terms and an afhne master equation [ ]. 

In the next section we will briefly discuss the pertur- 
bative open-system master equation which we will use to 
approach these issues. Then in Sees. Ill and IV we will 
provide resolutions based, respectively, on balancing the 
coupling and the correlations. 

II. OPEN-SYSTEM DYNAMICS 

In the time-local representation (also called the con- 
volutionless or Markovian representation), the dynamics 
of the reduced density matrix of the system p can be 
expressed with a quantum Liouville equation 

j^p{t)^C{t)p{t), (5) 

for any factorized initial condition. As a perturbative ap- 
proximation, C{t) is expanded in powers of the system- 
environment interaction Hi {t) and truncated to some or- 
der. (We momentarily consider full time dependence in 
this section as the more general relations will be neces- 
sary for our techniques.) Such perturbative master equa- 
tions can be derived in a variety of ways [ ~ ] and find 
application in many branches of physics and chemistry 
[10-13]. The expansion of C{t) will then take the form 

oo 

C{t) = Y,C2k{t), (6) 

*;=0 

£oP= HHs(<),p] , (7) 

where C2k — ©(Hj*^) and to zeroth-order the system is 
driven in a unitary manner by its Hamiltonian Hg [t). We 
take the expansion to be even as we only consider Gaus- 
sian noise, which is symmetric. A Gaussian noise distri- 
butional is necessary for higher-order perturbation the- 
ory to be non-secular in the late-time limit [ ] , though 
at second order one is effectively truncating the noise cu- 
mulants in a manner consistent with Gaussian noise. 

The linear master equation is derived under the general 
assumptions of a factorized initial state and an expansion 
of the interaction as a sum of separable operators: 

Hi(t)=^L„(t)®l„(t), (8) 

n 

where L„(t) and l„(t) are system and environment oper- 
ators respectively. The environment coupling operators 
\n{t) will typically be collective observables of the envi- 
ronment, with dependence upon very many modes. 

Using the notation of Ref. [14], the second-order master 
equation can be expressed as 

^2{P} = ^ [Ln, P (A„„0 L^)^ - (A„„0 lira) p] , 

nm 

(9) 



where the A operators and o product define the second- 
order operator 

(A„„oL„i)(t) = / dTanm{t,T) {Go(<,r)L„i(T)} . 
Jo 

(10) 

Here Gq (t,T) : p{t) — > p{t) is the free-system propaga- 
tor, which for a constant Hamiltonian H is given by 

Go(i, t)p= e-'{t-r)li p^+^(t-r)ll ^ 

while a„„i(i, t) are the environment correlation functions 
defined by 

a™(i,r) = (l„(t)Ur))E , (12) 

where l„(i) represents the time-evolving 1„ in the interac- 
tion (Dirac) picture. In general the correlation function 
is Hermitian and positive definite. For constant coupling 
to any stationary environment, the correlation function 
will also be stationary, q:(<,t) = a(t — T). Furthermore, 
for a thermal environment the correlation function will 
satisfy the KMS relation [L'j, IG]: 

a(cj) = a*(-cj)e-'^", (13) 

where aiuj) = dt a.{t) e^^'^'^ denotes the Fourier 

transform. 

From this perspective, the mathematical cause of the 
initial jolt becomes clear. For constant Hamiltonians and 
an initially stationary environment, the second-order op- 
erator obeys the relation 

^(A„„oL„)(i) = anmit) {Go(t)L„} , (14) 

which can be extremely large near the initial time when 
considering an environment with a sufficient amount of 
high frequency modes (such as low-temperature ohmic 
and supra-ohmic environments) since ck{t) is typically a 
very localized distribution in those cases. For a finite but 
large cutoff A, the correlation function becomes of order 
A for a time of order 1/A. 

III. COUPLING SWITCH-ON 

One method for balancing the initial coupling between 
the system and environment with their initial lack of cor- 
relation, is to turn on the coupling slowly with a time- 
dependent interaction such as 

Hi = 0,(f)^L„(g)l„, (15) 

n 

where 9s{t) : [0,oo) — > [0, 1) is a smooth switch-on func- 
tion with a characteristic timescale Tg, which vanishes at 
the initial time and becomes (effectively) one for times 
longer than Tg. To some extent, this was considered for 
linear systems in Ref. [ '"]■ 
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Such a time-dependent interaction is equivalent to em- 
ploying the second-order operator 

(A„„oL„)(i) = 6's(t) /dT6's(i-T)a„™(r){Go(T)L„} , 
Jo 

(16) 

for otherwise constant couplings and Hamiltonians. 
Therefore, any initial jolt due to the localized nature of 
a{t) will be suppressed by 9s{t) as long as Ts » 1/A. As 




nt 



FIG. 1. Zero-temperature, ohmic decay rate for the • instan- 
taneously coupled and ■ gradually coupled initial states of a 
two-level system with exponential cutoff frequency A = 100 SI. 
In this case the switch-on function is exponential, 9a{t) = 
1 — 6^'''^% and the switch-on times are chosen to take the 
values 1/A, 2/A, 4/A, 8/A, 16/A. 

can be seen in Fig. 1, the cutoff-frequency jolts are es- 
sentially replaced by jolts of frequency min[A, 1/ts] and 
amplitude proportional to the same value. This ap- 
proach provides a useful way of generating initial system- 
environment correlations when Tg is much larger than 1 /A 
but smaller than any other relevant timescales (such as 
the system frequencies). Furthermore, even if a mild jolt 
is still present, the important point is that it is cutoff 
insensitive (for fixed Tg and sufficiently large A). 



IV. DYNAMICALLY PREPARED INITIAL 
STATES 

Alternatively, in order to balance the initial correla- 
tions with an initially non- vanishing interaction strength, 
we will consider here initial states with suitable correla- 
tions to the environment. Such states will be obtained 
via an auxiliary construction which involves evolving an 
initially uncorrelated state for a sufRciently long time (a 
similar procedure was used in Ref. [18, 1!)] within the 
context of semiclassical gravity) . The system-enviroment 
correlations are then dynamically generated through the 
environmental interaction itself. Our first examples of 
equilibrium preparation will be the simplest mathemati- 
cally, while the final examples of non- equilibrium prepa- 
ration will be closer to actual laboratory experiments. 



In all cases we will take the system and environment 
to be uncorrelated not at t = but in the infinite past. 
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for some (possibly unimportant) system state pg{—oo) 
and thermal p-^{~-oo). We then define the system Hamil- 
tonian piecewise in time 



H m - ^ H+(t) < t 



(18) 



such that in past the system is allowed to equilibrate with 
the environment for an infinite time, which determines 
the correlated initial state at t = 0. The second-order 
master equation is then determined by 

(A„„iOL„)(t) = / dra„TO(t, r) {Gs(i, r)L™(T)} . 

J—oo 

(19) 

To analyze the coefficients associated with the initially- 
correlated state, we will reduce them to a sum of coeffi- 
cients for the auxiliary initially-uncorrelated state involv- 
ing various time ranges. First, we split the integration 
into two parts 



dT 



dr 



dr , 



(20) 



with the first integral depending only upon G+(t, r) and 
corresponding to the uncorrelated coefficients. Inserting 
the product G_ (0, t) G_ {t, 0), which equals the identity, 
the second integral can be written as 



M{t) / dr {t,T) {G_(i,T)L„(T)} , 

oo 



(21) 



given the operator 

M{t) = G+{t,0)G^{0,t). (22) 
The integral in Eq. (21) is then broken up into two parts 



pt 

dr = dr 

OO J—OQ 



dr , 



(23) 



corresponding to the asymptotic and finite-time coeffi- 
cients for an initially uncorrelated system driven by the 
time-independent preparation Hamiltonian H_ . Finally, 
our correlated coefficients can be expressed in terms of 
the uncorrelated coefficients as 



(A„„oL„)(t) = (A„,„oL„) , (i) 



(24) 



correlated 



uncorrelated 



Al(t)|(A„„oL„)_(i) - (A„„oL„)_(oo)| 



jolt suppression 



preparation eraser 



where the subscripted (AoL)_|_ coefficients are defined 



as 



(A„„iOL„)±(t) = / dTa„m{t,T) {G±(t, T)L,„(r)} 

"'0 



(25) 
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If the system frequencies are always small as compared to 
the cutoff, we can inspect the early-time behavior (and 
jolts) by letting G±{t) « 1. Then one can see that the 
first two terms of Eq. (24) will precisely cancel in the 
early-time regime. Therefore, the correlated initial states 
are jolt-free given sufficiently small system frequencies as 
compared to the cutoff: 17 ^ A. The final term in turn 
is such that in the late-time limit it precisely cancels the 
second term and erases all memory of H_ . Finally note 
that, quite trivially, if we choose H_|_(i) = H_, then the 
first two terms cancel and we recover the equilibrium 
coefficients at any finite time. 



A. Equilibrium preparation 

To prepare an initial state in this approach, we choose 
the past Hamiltonian H_ such that its dynamics along 
with the environment interaction relaxes our system to 
the desired initial state: 

lim e*^-(°°)po-Po, (26) 

r— >-oo 

£_(c^)Po = 0, (27) 

where C-{oo) is the stationary limit of the Liouvillian 
for a system with the past Hamiltonian as well as the 
coupling to the environment. 

Our target state Pq will only be specified to zeroth 
order in the system-environment interaction. This is be- 
cause for sufficiently long times (and in particular for the 
asymptotic equilibrium state) the diagonal elements of 
the reduced density matrix in the energy basis cannot be 
determined beyond zeroth order anyway when using the 
second-order perturbative master equation [20]. Due to 
unavoidable degeneracy present in all open-system dy- 
namics, one actually requires components of the fourth- 
order master equation to calculate the full second-order 
solutions. The second-order master equation provides for 
all second-order dynamical quantities, such as frequency 
shifts, dissipation, diffusion and decoherence rates. We 
are concerned here with the induced jolts, which are dy- 
namical quantities, and so this subtle point does not raise 
any additional problems for us. 



1. Preparation by decoherence 

For L„ all commuting with each other, one can force 
a general environment into ^-state preparation via de- 
coherence. If the past Hamiltonian is deactivated, or 
more generally taken to commute with L„, then since all 
system operators commute with each other, the master 
equation and its solutions will trivially result in a system 
which decoheres in the ^-basis associated with the L„. 
Thus, coefhcients prepared in this manner are consistent 
with any initial state which is a completely incoherent 
mixture of ^-states. [Note that if pg(— cxo) corresponds 



to a pure eigenstate of the set {L„}, this procedure sim- 
ply adjusts the state of the environment, while system 
and environment remain unentangled.] 



2. Preparation by equilibration 

A finite-temperature environment allows mixed state 
preparation by equilibration. Essentially one chooses the 
past Hamiltonian so that its thermal state (or some other 
steady state) is the desired initial state. For a positive- 
temperature environment, at zeroth order one can pre- 
pare a (sufficiently) mixed state Pq with the past Hamil- 
tonian H_ = — T log(pg). However, one must be careful 
that past system frequencies are small as compared to the 
high frequency jolts, otherwise this preparation will fail 
to remedy jolting. One can work out that the adiabatic 
preparation regime is given by 

^«e^^, (28) 

Pmin 

where A is the jolt frequency and p are the initial state 
probabilities of preparation energy levels connected by 
L„. (Clearly, for this method to work there can only be 
a finite number of such energy levels.) 



3. Preparation by freezing 

To prepare an initially pure state via equilibration 
at the order that we are working, one requires a zero- 
temperature environment for preparation by freezing. 
Then one can choose any H_ with ground state pg. It is 
important to emphasize that the reduced density matrix 
of the system corresponding to the ground state of the 
combined system will not be a pure state in general due 
to the entanglement between the system and the environ- 
ment: the free ground state of the system is a pure state, 
but the reduced density matrix of the open system is in 
general a mixed state beyond zeroth order in the system- 
environment coupling. However, this point becomes irrel- 
evant at the order that we are working since, as explained 
above, when using the second-order perturbative master 
equation to prepare the initial state by equilibration, one 
cannot meaningfully specify Pq beyond zeroth order. 



B. Non-equilibrium preparation 

In order to consider situations closer to actual lab- 
oratory experiments, here we will first allow the sys- 
tem to equilibrate with the environment (as described 
in the previous subsection) and then choose some prepa- 
ration Hamiltonian Hp(i), which would (in the absence 
of coupling to the environment) generate the desired ini- 
tial state in some finite time rp. One simply applies 
the master-equation coefficients in Eq. (24) with future 
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while evolving the combined system during this addi- 
tional finite preparation time. In fact, the preparation 
time Tp cannot be too long if we want the state of the 
system to be more or less close to lipo). 



2. State swapping 



FIG. 2. Zero-temperature, ohmic decay rate for the • un- 
prepared and • prepared initial states of a two-level system 
with exponential cutoff frequency A = 100 Q. In this case 
preparation by freezing was used to create an initially excited 
state. 



Hamiltonian 



Ho(i) 
Hp(t) 



Tp < t 
t < Tp 



(29) 



where Ho(i) is the desired post-preparation Hamiltonian. 
All jolts will be avoided if 1/rp ^ A: the introduction 
of a non-vanishing preparation frequency serves to tame 
the jolts and eliminate their high-cutoff sensitivity. 



Let us consider a system which is initially equilibrated, 
without making any assumption as to the temperature of 
the environment. We couple the system to an ancillary 
and analog system (equivalent Hilbert spaces) that is al- 
ready prepared in the desired initial state. The system 
of interest and ancilla are temporarily coupled in such a 
way that they swap states, for instance by means of the 
following block-matrix preparation Hamiltonian: 



Hp 



TT 

2t> 



1 

1 



(31) 



In the absence of coupling to the environment this would 
exactly swap the system and ancilla states in a time Tp. 
The same remarks as for state flipping concerning the 
purity and accuracy of the prepared state when taking 
into account the coupling to the environment also apply 
in this case. 



3. Other possibilities 



1. State flipping 

A possible preparation Hamiltonian, which could 
model as a particular case Rabbi oscillations induced by 
an appropriate laser field acting on a two-level system, is 
the following: 

Hp = ^ (IV'o)(0| + lOXV^ol) . (30) 

Assuming that one has a zero-temperature environment 
and that the system is already equilibrated, driving the 
system with this Hamiltonian for a time Tp provides a rel- 
atively easy way of preparing an initial pure state 
As discussed above, the reduced density matrix of the 
system will actually be a mixed state in general, because 
of the system-environment entanglement of the equilib- 
rium state as well as the interaction to the environment 



Within the second-order perturbative approach, gener- 
ation of equilibrium correlations in a laboratory setting 
can always be calculated using Eq. (24). One only needs 
to make sure that any additional state preparation does 
not rely upon large system energies as compared to the 
bath cutoff. For instance, one can consider the prepa- 
ration of Ref. [5], which relies on ancillary degrees of 
freedom to drive the equilibrium state into a coherent 
superposition. In fact, one could simply apply their own 
time-dependent Hamiltonian to our formulas as H+(t) 
and obtain results consistent with theirs. 
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